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The properties of the QCD partition function at finite chemical potential are studied 
within the instanton liquid model. It is shown that the density dependence of the quark- 
^ , induced instanton-antiinstanton (I-A) interaction leads to the formation of topologically 
\^ [ neutral I-A pairs ('molecules'), resulting in a first order chiral phase transition at a critical 



chemical potential jj^ ~ 310 MeV. At somewhat higher densities (/Xg > 360 MeV), the 

quark Fermi surface becomes instable with respect to diquark condensation (Cooper pairs) 

O . generating BCS-type energy gaps of order 50 MeV. 
OO 

a^ 

1. INTRODUCTION 

The investigation of the phase diagram of QCD is one of the central issues in under- 
oy, standing the properties of strong interactions. While the finite temperature axis has been 



(N 



^ ' theoretically explored in quite some detail, much less is known about the finite density 
(/ig-) axis, which is partly due to the fact that first principle QCD lattice calculations at 
^ ■ finite fig encounter the problem of a complex fermionic determinant when integrating the 
QCD partition function. The phase structure of QCD at finite density, however, might 
be very rich, including new forms of condensates (other than the ordinary chiral conden- 
sate present in the QCD vacuum), different orders of phase transitions in the /Xg-T-plane 
(entailing a tricritical point |jl|), etc.. In this contribution, however, we will focus on the 
T=0-, fiq>0- axis and try to examine the nature of the chiral phase transition in this 
regime. We will do so within the framework of the instanton liquid model (ILM) of QCD, 
which, although lacking explicit confinement, yields a very successful phenomenology of 
the QCD vacuum structure, and the low-lying hadron spectrum. Recently it also provided 
an interesting mechanism for chiral symmetry restoration at finite temperature, based on 
the rearrangement of the (anti-) instantons within the liquid, rather than on a mere dis- 
appearance of them as suggested earlier. Our objective here is to investigate whether a 
similar mechanism could be responsible for the restoration of chiral symmetry at finite 
density as well. 

We start by briefly recalling some features of the instanton liquid model at zero density 
(sect. 2) and then turn to the finite density case (sect. 3), where we first calculate the /in- 
dependence of the instanton-antiinstanton interaction and then assess its impact on the 
chiral phase transition. Under certain approximations we will be able to avoid a complex 
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partition function and estimate the critical chemical potential within a mean-field type 
approach. We furthermore include some correlations between quarks at the Fermi surface 
and show how the recently discussed mechanism of color superconductivity figures 
in our mean-field description. 

2. THE INSTANTON LIQUID MODEL AT ZERO DENSITY 

The instanton model in QCD is based on the assumption that the vacuum is dominated 
by non-peturbative gauge field configurations which constitute semiclassical ('instanton') 
solutions of the Yang-Mills equations. The QCD partition function in the instanton ap- 
proximation then becomes 

Zqcd= E n dn, nipj) e~'"- pf^ U det{i P + tnif) , (1) 

N+,N- ^^ + -^^~- 1=1 f=i 

where the path-integral over all possible field configurations has been converted into an 
integration over the so-called collective coordinates flj = {zi,pj,uj} (position, size and 
color orientation) of iV_|_ instantons and N_ antiinstantons. The single-instanton am- 
plitude n{pi) (including quantum corrections) and the gluonic part of the instanton in- 
teraction Sint determine the total instanton density in the absence of ferminons to be 
N/V = 2 J dp n{p) e~'^'"*. The determinant of the Dirac operator, arising from the inte- 
gration over the quark fields, is approximately calculated by keeping only the lowest lying 
modes, which should be reasonbale for assessing the bulk properties of the system. The 
Dirac equation in the (anti-) instanton field possesses a (right-) left-handed zero energy 
solution. In the basis spanned by these zero modes, \E'o,/ and ^E'o.a, and neglecting small 
current quark masses {rrif 0), the fermionic determinant reads 



t/et(^ ^) ~ rfet J^j=|T,Ar, (2) 

where the overlap matrix element 

Tja{z, u) = J d^x '^lj{x - zi, ui) {i JP) '^\^a{^ - za,ua) = i u- z f{z) (3) 

is linear in the relative S'?7(3)-color orientation characterized by a complex four vector 
{z = za — zi denotes the relative distance between J and A). Lorentz invariance then 
implies that Tia is determined by a single scalar function /(v^)- 

With the key parameters taken as N jV ~ (1-1.4) fm~^ and pi = pA— 1/3 fm (which 
have also been confirmed in lattice calculations) a successful phenomenology of the QCD 
vacuum and the low-lying hadronic spectrum can be obtained [Q. 

3. THE INSTANTON LIQUID MODEL AT NON-ZERO DENSITY 

3.1. 1-A Interaction at Finite pq. Quark Zero Modes and Tia 

To study medium modifications of the instanton ensemble we first have to construct 
their interactions at finite quark chemical potential. We start from the finite density Dirac 
equation (in euclidean space), which still has zero mode solutions satisfying 

(z^j-i/ig74)^o,/ = . (4) 
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Constraining ourselves to zero temperature, the gluonic (instanton-) fields entering the 
covariant derivative are not affected by fig. The explicit form of the quark zero modes has 
been determined in ref. 0, 



with the function Il{x) = (1 + p'^/x^) and a spin-color spinor x- Note that the solution of 
the adjoint Dirac equation, 

*o,/(a;; -f^q) - *^g74) = , (6) 

carries the chemical potential argument with opposite sign. This is necessary for a con- 
sistent definition of expectation values at finite fig, and in particular renders a finite norm 

J d^x ¥q/x; -fig) ^o,/(a;; Aig) = 1 . (7) 

With the properly constructed quark wave functions we can now evaluate the fermionic 
overlap matrix element representing the zero mode part of the full Dirac operator. Choos- 
ing for simplicity the sum ansatz for the gauge-field configurations, = Ajj^+A^, entering 
the covariant derivative, allows us to replace the latter by an ordinary one, yielding 

Tja{z, u; fig) = -J d^x '^lj{x - zi; -fig) (i ^ - ifigli) ^0,^(2; - Za\ fig) 

= i U4 fi{T,r;fig) +i —— f2{T,r;fig) . (8) 

The breaking of Lorentz invariance in the medium implies the existence of two independent 
(real- valued) functions /i,/2 which are shown in fig. 1. Similar to what has been found 
at finite temperature [^], Tja is strongly enhanced in temporal direction with increasing 
fig, but oscillates as ~ sin(2/iqr) in spatial direction 0. The latter effectively suppresses 
the interaction once integrating over r. 

3.2. Thermodynamics of the ILM at Finite fig 

To investigate the finite density properties of the thermodynamic potential (or free 
energy) we here resort to the so-called cocktail model introduced in ref. . It amounts to 
a mean-field type description including three major components in the system: essentially 
random (anti-) instantons (the 'atomic' component), strongly correlated /-A-pairs (the 
'molecular' component) and a Fermi sphere of constituent quarks ('quasiparticles'). Thus 

n{flg) = Qinstifig) + ^guarkiN) ' (9) 

where the constituent quark contribution 

^quarki^q', = ^qi^q'^ f^q) " f^q ngivUg, fig) (10) 

is zero for Fermi energies fig < irig, with irig denoting the constituent quark mass. The 
instanton part of the free energy. 
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Figure 1. Quark- induced /-A- interaction at finite density for tlie most attractive color 
orientation U4=l, u=0 as well as r=0 (left panel) and for ^4=0, \u\=l and r=0 (right 
panel) . 



is related to the atomic and molecular 'activities' 

= p2(^-^) e-2^- {[TMr^Mr'^f) . (12) 

The minimization of Q with respect to the corresponding (4-dimensional) densities Ua 
and rim determines the equilibrium state of the system at fixed /ig. In particular, Qmst 
encodes all the features of the T = fig = instanton ensemble, as e.g. the quark con- 
densate and the constituent quark mass, which in mean-field approximation are given 
by {qq) = — l/(7rp)(3/2na)^/^ and nig cx —p'^{qq). Therefore the normalization constant 
C oc {AqcdY can be fixed to give N/V=1.4 fm~^, being realized for na=1.34 fm~^ and 
nm=0.03 fm~^, which is not unreasonable. The gluonic interaction has been approximated 
by an average repulsion Sint = —i^p'^{na + 2nm). 

At finite pg the Dirac operator is not hermitian any more, i.e. TAi{pq) 7^ Tj^{pg), re- 
sulting, in general, in a complex fermionic determinant, entailing the well-known 'sign' 
problem. However, assuming an average gluonic interaction that does not depend on 
density allows us to perform the color averages implied in eqs. (0) analytically, i.e. 



Za OC J duTiA{pq)Tj^{-pg) = ^[f+f\ + f+f^ 



Zm cc J du[TiA{Pq)TlA{-pg)] f = 4NciN^ - 1) ^ ' 

{Nf = 2, = /i(±yUg)), which ensures the pressure to remain real. Fig. 2 shows our 
results as function of the quark chemical potential (left panel) for Nc = 3, Nj = 2. At 
small pq essentially nothing happens until, at a critical value /i^ ~ 310 MeV, the system 
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Figure 2. Our results for the instanton ('atomic') and molecule densities (upper left panel), 
the pressure p = —Q (middle left panel) and the constituent quark mass (lower left panel) 
after minimizing the free energy, eq. @, with respect to Ua and rim] additionally including 
scalar diquark correlations (as discussed in sect. 3.3) adds the dashed-dotted line to the 
lower left panel, representing the BCS energy gap at the quark Fermi surface. The right 
panel shows the free energy as a function of constituent quark mass niq oc nl^"^, indicating 
a first order transition from the minimum at finite rriq to the one at rriq = 0. 



jumps into the chirally restored phase, the latter being characterized by = 0. The 
transition is of first order, as can be seen by inspection of the mg-dependence of the 
free energy (right panel of fig. 2). Below /i^, the pressure actually decreases slightly with 
increasing fiq indicating a mixed phase-type instability, similar to what has been discussed 
in refs. @,|^. A significant difference, however, is given by the fact that in our approach 
the total instanton density at the transition (residing in J-A-molecules) is still appreciable, 
N/V = 2nm — 1.1 fm~^, providing the major part of the pressure at this point; in other 
words: a substantial part of the nonperturbative vacuum pressure persists in the chirally 
restored phase (of course, eventually it will be suppressed due to the Debye screening of 
the instanton fields, which we have not included here). 
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3.3. Color Superconductivity 

As has recently been pointed out in refs. 0,^, the quark Fermi surface in the plasma 
phase might be unstable with respect to the formation of quark-Cooper pairs, once an 
attractive q-q interaction is present. Using the instanton- induced interaction in the scalar, 
color-antitriplett diquark channel (which is essentially the Fierzed-transformed interaction 
leading to a deeply bound pion and is also phenomenolgically well-supported by baryonic 
spectroscopy), BCS-type energy gaps of Aq — 50-100 MeV have been predicted. 
In the mean-field approach employed here, this interaction modifies the quark quasipar- 
ticle contribution according to 

= tr log [DK, Ao)] - tr [Z^K, Ao)S(Ao)] + Gtr [FK, Aq)] tr [fK, Aq)] (14) 

with the quasiparticle quark-propagator D and the annomalous (Gorkov) propagator F 
corresponding to creation/annihilation of a bound Cooper pair. G is the effective coupling 
constant derived from the instanton g-g- vertex . From one can obtain the standard 
gap equation for Aq. However, here we also have to account for the fact that a finite 
Aq will damp the quark zero-mode propagators in the instanton part of the free energy, 
which results in a suppression of the quark induced I- A interaction Tja of eq. (H). Thus, 
^inst disfavors finite Aq's. The resulting expression for the free energy is of the form 

Q{flq) = Qinst{na,nm,Ao; fig) + ^ 2Qf {iTlq, Aq] fig) + {N^ - 2)Qq{mq; fig) , (15) 

(the last term accounting for unpaired quarks), which now has to be minimized w.r.t. Ua, 
rim and Aq. We find that color superconductivity does not appear before fig ~ 360 MeV 
(see dashed-dotted curve in the lower left panel of fig. 2): although the quark-part of the 
free energy by itself, eq. (p!4D, always favors a finite value for Aq, the suppression caused 
by Aq in flmst {i-e. the damping of quark-propagation in Tja) prevents the formation of a 
diquark condensate below fig ~ 360 MeV. This again reflects the fact that the instanton 
contribution to the pressure is still dominant in the region somewhat above /x^, so that the 
gain due to a finite Aq in fl'^ cannot overcome the 'penalty' in flinst- Above fig ~ 360 MeV 
energy gaps of order ~ 50 MeV arise, in line with the findings of ref. 0. 
As far as the thermodynamic properties of the system are concerned, no changes as 
compared to the previous section occur below fig = 360 MeV, and even above the results 
for na{fig), rimifiq), pifJ'g) and nig^fig) would be hardly distinguishable from the curves 
displayed in fig. 2. 



In summary, we have studied the QCD partition function at finite density using the 
instanton model. Employing a simple mean-field approach, chiral symmetry restoration 
emerges at fig =310 MeV as a first order transition from an essentially random (anti-/) 
instanton liquid into a phase of strongly correlated I- A molecules (and massless quarks), 
driven by the density dependent increase of the quark-induced /-^-interaction. Neglecting 
possible medium modifications in the gluonic interaction, which are expected to be small, 
a complex pressure could be avoided. Accounting for correlations in the scalar diquark 
channel, color superconductivity sets in for chemical potentials fig > 360 MeV, associated 
with BCS gaps of up to ~ 50 MeV. 
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